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RECIPROCITY RELATIONS IN AERODYNAMICS! 


By Max. A. HEAsLET and Joun R. SPREITER 


SUMMARY 


Reverse-flow theorems in aerodynamics are shown to be based 
on the same general concepts involved in many reciprocity the- 
orems in the physical sciences. Reciprocal theorems for both 
steady and unsteady motion are found as a logical consequence 
of this approach. No restrictions on wing plan form or flight 
Mach number are made beyond those required in linearized 
compressible-flow analysis. A number of examples are listed, 
including general integral theorems for lifting, rolling, and 
pitching wings and for wings in nonuniform downwash fields. 
Correspondence is also established between the buildup of circu- 
lation with time of a wing starting impulsively from rest and 
the buildup of lift of the same wing moving in the reverse 
direction into a sharp-edged gust. 


INTRODUCTION 


Some of the most important results in the recent study 
of wing theory have been achieved through the development 
of reverse-flow relations. The theorems already obtained 
are of outstanding practical utility and it appears obvious 
that the fullest exploitation of the methods has yet to be 
accomplished, either from a purely theoretical standpoint 
or in the routine calculation of wing characteristics. Atten- 
tion. to such problems in aerodynamics was initiated by von 
Kármán (ref. 1) who first announced the invariance of drag 
with forward and reversed directions of flight for a nonlifting 
symmetrical wing at supersonic speed. Subsequently, 
advances in the theory were made by Munk, Hayes, Brown, 
Harmon, and Flax (refs. 2 through 7). Up to the present 
time, the most general results have been expressed by 
Ursell and Ward (ref. 8) and by R. T. Jones (ref. 9) in his 
attack on the study of wing shapes of minimum drag. 
In the forms given in the two latter papers, the derived 
equivalences could be termed reciprocal or reciprocity rela- 


tions rather than reverse-flow relations; in fact, this change 


in terminology divorces attention, momentarily, from the 
purely aerodynamic aspects of the results and, in this way, 
suggests reorientation in terms of the various similar 
relations appearing in other engineering fields. In the 
theory of elasticity, for example, a reciprocity theorem for 
small displacements of an elastic medium is so expressed 
as to appear in formal agreement with the statement of the 
result given by Ursell and Ward (see, e. g., ref. 10). This 
theorem is attributed to E. Betti and was published in 1872. 
A generalization was given by Lord Rayleigh in 1873, and 


in various sections of his two volumes on the theory of sound 
(ref. 11) discussions of ‘reciprocal relations in an elastic 
medium and for acoustic sources are given. In 1886, von 
Helmholtz (ref. 12) obtained, by means of variational 
methods applied to Hamilton’s characteristic function, a 
reciprocal theorem for small changes in the momenta and 
coordinates of a general dynamical system in forward and 
reversed motion. This result was commented on, in turn, 
by Lamb (ref. 13) and an independent proof based upon 
Lagrange’s generalized equations was given. The paper by 
Lamb is of particular interest since it contains the essential 
idea underlying the development of reverse-flow theorems 
in wing theory. Thus, Lamb remarks, as had Lord Rayleigh 
previously, that reciprocity relations between sound sources 
do not apply directly in a moving atmosphere. He points 
out, however, that the reciprocity can be restored if the 
direction of the wind is also reversed. Further examples 
of reciprocal theorems appear in the theories of electricity 
and magnetism (in particular, reference should be made 
to Maxwell’s discussion of the subject in ref. 14) and of optics. 

The generality in the statement of reciprocity relations 
appears, almost universally, to have held back their applica- 
tion to problems for which they are obviously, in retrospect, 
particularly fitted. This generality is even more apparent 
in some of the conclusions of Lord Rayleigh and von Helm- 
holtz which apply to nonconservative systems. 

The purpose of the present paper is twofold. First, a close 
connection. will be established between reverse-flow theorems 
in subsonic and supersonic, steady-state wing theory and 
known reciprocity relations between two solutions of the 
equation governing the flow field. In this way, machinery 
will be provided whereby extensions of existing results to the 
case of unsteady motion follow directly. Second, a number 
of particular problems in wing theory in steady and unsteady 
flow will be considered. It will be shown that, provided 
attention is limited to force and moment characteristics, 
the complexity of many solutions involving nonuniform: 
flow fields, control-surface deflections, and unsteady motion 
can be reduced considerably. In some cases, previously 
obtained solutions will be calculated. Comparison with 
the original calculations will almost invariably highlight 
the economy of. effort in obtaining the final result. The 
utility of reverse-flow theorems is based on the fact that 
they build from known solutions and thus avoid the necessity 
of starting each problem anew. 


I! Supersedes NACA TN 2700, “Reciprocity Relations in Aerodynamics” by Max. A. Heaslet and John R. Spreiter, 1052. 
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GENERAL ANALYSIS 


RECIPROCITY RELATIONS FOR A CLASS OF PARTIAL 
DIFFERENTIAL EQUATIONS 


. In this section, integral relations associated with linear 
partial differential equations will be reviewed from the stand- 
point of relating independent solutions. The subject 
matter is precisely Green’s theorem and, in common with 
the usual expression of the theorem, it is preferable to 
treat the variables initially as abstract quantities. Consider, 
therefore, a class of linear partial differential equations of 
second order with independent variables Xi, Xa, ..., Xm 
that may be thought of as rectangular coordinates in a 
space of m dimensions. Denote differentiation of the 
function Y(X, Xz. . - , Am) with respect to the variables 
Xx; and X; by the subscript notation 


oy 
(Y)y= OX10.4; 


and consider the differential equation 
Ly)=}; ey A) + By=0 (1) 
i=] j=l 


where, for the purposes at hand, Ay,y—A;; are made inde- 
pendent of X, and X;, and Bis a constant. Such equations 
fall within the class of self-adjoint equations. 

By Green’s theorem (see, e. g., ref. 15), it is possible to 
relate two arbitrary functions y and © by the integral expres- 
sion 


ff form-o = f f D, y— yD, dS (2) 


where the left member is a volume integral over a prescribed 
region in m-dimensional space and the right member is a 
surface integral over the hypersurface S enclosing the given 
region. Equation (2) certainly holds for any region in 
which y and Q and their first and second derivatives are 
continuous. The directional derivaties D, are defined in 
terms of the direction cosines 7, na, .. -, Mm of the normal to 
the surface S with the stipulation that the normal is directed 
into the given region. Thus, setting 


m a 


2) My No 


where vı, ¥2,.--, Vm are the direction cosines of a line termed 
the ‘‘conormal,’” the directional derivative is defined by the 
expression 


Div=N SNN =N O 


If, finally, y and Q are assumed to satisfy equation (1), the 
left side of equation (2) vanishes and the resulting expression 


f f OD, pd S= f f YD.2AdS - (4) 


is a general reciprocity relation expressing the functional 
dependence between two arbitrary solutions of equation (1). 

An interesting interpretation of equation (4) has been 
given by several writers (see, e. g., ref. 16, p. 46) and applies 
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to the particular case when y and 2 are identified with the 
perturbation velocity potential » in the theory of incom- 
pressible-fluid flow. The governing equation is Laplace’s 
equation in three dimensions 


Pret Pyr F Yr2=0 (5) 


and the reciprocity relation takes the form, for two possible 
solutions ¢ and g’, 


ffe £ as=f fe $2 as (6) 


where the directional derivatives are now along the true 
normals to the surface © enclosing a three-dimensional 


volume. It is known that any actual state of motion of a 
liquid for which a single-valued velocity potential exists can 


be produced instantaneously by the application of a properly 
chosen system of impulsive pressures. These impulsive 
pressures are directly proportional to the velocity potential 
plus an arbitrary constant which: may, in the present case, 
be associated with the pressure of the uniform stream. 
Equation (6) is thus seen to represent summations of the 
cross products of impulse and normal velocity in two possible 
motions of a conservative system and is a special case of the 
dynamical theorem (ref. 11, p. 98) 


> Prr = > Pr Qs (7) 
Teal Tml 


where p,, g, and p,’, g,’ are generalized components of 
impulse and velocity in any two possible motions of a system 
which starts from rest. 

The interpretation of equation (6) that leads to equation 
(7) provides an indication of the close connection between 
reciprocal theorems based upon the principles of least action 
and the symmetric character of Green’s theorem for certain 
second-order differential equations. In the subsequent 
applications it will be convenient to proceed directly from 
equation (4) and seek to establish reciprocal relations between 
flow fields in wing theory. Such a process is well known 
when y or Q is replaced by the elementary solution associated 
with a unit source and, in this case, establishes a general 
solution in terms of source and doublet distributions deter- 
mined by arbitrary boundary conditions.: The present 
objective is, however, different in that one wishes to get a 
symmetrical dependence between two general solutions. 
Moreover, the apparent symmetry of equation (4) must be 
consistent with physical considerations so that, for example, 
the velocity at point A induced by a single source at point B 
is equal to the velocity at point B induced by a source at 
point A. If, in the two systems, the effective free streams or 
flight directions are opposed, a fore and aft symmetry occurs 
and the possibility of maintaining symmetry in the reciprocal 
relations becomes feasible. 


RECIPROCAL THEOREMS IN WING THEORY 
REVERSE FLOW FOR SUBSONIC WINGS 
Consider a thin wing at possibly a small angle of attack and 
situated in the immediate vicinity of a flat surface Which is 
designated the zy plane. For sufficiently small thickness and 
angle of attack, the perturbation velocity potential or any of 
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the perturbation velocity components of the wing satisfy the 
linearized partial differential equation of compressible flow. 


Bee Vay t Yrr= 0 (8) 


where 6?=1—M,?=1—(Up/ao)?, and Uco are, respectively, 
the flow velocity and speed of sound in the free stream. 
Equation (8) applies in forward or reverse flow, provided the 
corresponding free-stream Mach numbers M, and My’ are 
equal, In figure (1) a lifting wing with plan form P is 


indicated along with the vortex wakes as the configuration 
would appear if the two flow fields were superimposed. It 
will be assumed that the wing chord is finite and that the 
profiles are closed. 


FIGURE. 1.-—-View of wing in combined flow field. 


It is proposed to apply the general reciprocity relation (4) 
to these flow fields in a manner similar to the approach used 
in developing the basic solutions of the differential equation 
(sco, e. g., refs. 17 and 18). Thus, for subsonic flow, hemi- 
spherical regions of large radius and lying first above and 
then below the plane of the wing will be chosen as the 
volumes of integration. The surface integrals will therefore 
extend over a hemisphere with center at the wing and a flat 
surface that lies immediately adjacent to the z=0 plane. 
This latter surface is subsequently to be brought into coin- 
cidence with the plane of the wing but must be considered 
first in its displaced position since only then can the flow be 
assumed free of possible singularities in perturbation veloci- 
ties and their gradients. As in wing theory, in general, the 
attenuation of the perturbation potential and its gradients 
may be assumed of such a nature that the integrated con- 
tributions of the wing and its wake over the hemispherical 
surfaces vanish in the limit as the radius becomes infinitely 
large. It remains, therefore, to consider the integrals over 
the surfaces at z=0-+e and z=0—e. Denote these surfaces 
as o, and o1, respectively, where the subscripts u and 1 specify 
values above and below the z=0 plane. Equation (4) can 
then be written as 


ff y SE dS = f f yp’ ov dS; -n directed upward (9a) 
O u - 


Ty 
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and 
ff Y ov dS = f y oe as; n directed downward (9b) 
O? Or 


where the primes denote conditions for the reversed flow. 

In equations (9), is now replaced by the perturbation 
velocity potential y(z,y,z) and ¥ by the z or streamwise 
component of perturbation velocity u(z,y,z). By virtue of 
the wrotationality of the flow, the gradients of u and w are 
related by the expression 


Ou _ 


dw 
n E 


Ox 


and equations (9) can both be written in the form 


[fo asf fees 


Integration by parts over either of the surfaces o, or o; gives 


f gw’ g -dy— fuwas=f fwwas 


At z=—o the potential ¢ for the forward flow vanishes and 
at z=-+o the upwash w’ in the reversed flow vanishes so 


(10) 


~ that the first term on the left is zero. , The remaining double 


integrals have for their surfaces of integration the displaced 
planes o, and o; In order to obtain a concise form of the 
reverse-flow theorem, it suffices to subtract the integrals ex- 
tended over. o; from the integrals over o„ and let the planes 
approach coincidence with the plane of the wing. Since 
ww and u,v’ are continuous everywhere except possibly 
in the immediate vicinity of the wing, the integration region 
can be restricted to planes slightly above and below the wing 
but extending beyond the wing edges. Provided the singu- 
larities at the edges can be disregarded, the analytical ex- 
pression of the steady-state reverse-flow theorem of Ursell 
and Ward (ref. 8) becomes 


-| { (u,w,’ —upo/)dS= f (Us Wy—U,/w)dS . (11) 
P P 


for either lifting surfaces or symmetrical wings where P is the 
plan form of the ‘ing in the z=0 plane. 

It remains to discuss the effect. of edge singularities. In 
the case of a lifting surface, square-root singularities in both 
u and w can occur at the leading and trailing edges just on 
and off the wing, respectively. In the combined flow fields, 
the limiting process would then yield residue terms analogous 
to the leading-edge thrust of a lifting -plate. If, however, 


_ the Kutta condition is imposed at the trailing edge for both 


the flow in forward and reverse direction, a combination of 
singularities does not occur and equation (11) is valid. If 
the leading edge of a symmetrical wing is rounded so as to 
produce a square-root singularity in w on the wing, a square- 
root singularity in u occurs just off the wing and a term cor- 
responding to the leading-edge drag (ref. 19) appears. If 
the geometry of the wing is fixed in the forward and reverse 
flow, however, the effect of these terms is canceled and equa- 


r 
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tion (11) still applies. It is important to point out that if 
the Kutta condition is not satisfied, then the area of integra- 
tion of equation (11) cannot be restricted to the plan form. 

The two sides of equation (11) are expressed in terms of the 
_ same coordinate system but it is usually preferable to asso- 

ciate with each of the two streams an x axis extending in the 
stream direction. To this end introduce now the subscripts 
1 and 2 to denote the forward and reverse flow and the two 
coordinate systems. Thus, in general, 


et (12) 


where ¢ and 7 are arbitrary constants, and equation (11) then 
becomes 


X= — r+ é, Yi=—Yrtn, 


f te (EY) Way EN) — Ur, CY), ry) andy = 


1 


f f Diu, (2,2) Wu, (ajfe) Un (asa) Wr, (ayy) Idexdy, (13) 


2 
In the case of a symmetrical nonlifting wing, the relations 


Uy, SU Wu = — Wy; 


must apply and, in the case of a lifting surface, linearized 
theory yields 


Uy™=— Un Wa =W] 


It follows that in eaa case, equation (13) reduces to the 
form 


ff UEY) Ww 1, yi) andy; = f fu (a2,Ya) Wi (£2,Y2)dzzdy2 (14) 
Pı P: . ` 


where the velocity components can be evaluated on either 
the upper or lower surface of the plan form. MH, further- 
more, the linearized pressure relation 


P—Do=— pou ' (15) 


is used where p is local static pressure, Po is static pressure 
in the free stream, and the wing profiles are assumed closed, 
equation (14) becomes i 


[J rrceadratensrddrdr=[ f penyo eydda (16) 
Py P: 


If, instead of specifying boundary conditions in a single 
plane, it is necessary to treat boundary conditions for a 
system of planes, the expression of the reverse-flow theorem 
is of the same general form as equation (11). Provided the 
Kutta condition is imposed at the trailing edges of all lifting 
surfaces, the relation becomes 


fS Poruds= ff Vaas (17) 
p> p>; 


where the area of integration 2 extends over both sides of 
all the wing surfaces, V, is the component of perturbation 
velocity normal to and directed away from each wing, and 
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the subscripts 1 and 2 refer to forward and reverse flow in 
the two axial systems. 


, REVERSE FLOW FOR SUPERSONIC WINGS 


The development of a reverse-flow theorem for supersonic 
wings parallels closely the analysis for the subsonic case. 
For either planar or multiplanar systems, the conormal in 
equation (4) is, in fact, the normal so long as tho surface of 
integration is a plane parallel to the v axis. In the case of 
the single wing, for example, equations (9) apply where the 
surfaces o are slightly removed from the plane of the wing 
and where y satisfies the differential equation 


(Mi— 1) Vor— Poy — Yee =0 


In the limit as e approaches the z=0 plane, the reversibility 
theorem takes the form of equations (11) and (16), provided 
the integration extends beyond the edges of the wing. It is 
necessary to include these edges for wings with subsonic 
leading and trailing edges since singularities occur in the 
perturbation components and.the solutions are not neces- ’ 
sarily unique. For supersonic-type edges, the area of 
integration can be confined to the plan form of the wing and 
this is also true for subsonic edges, provided the Kutta con- 
dition holds for all subsonic trailing edges in both the forward 
and reverse flow. Equation (17) relates the two possible 
flows in the case of multiplanar systems. 


REVERSE UNSTEADY MOTION 


In the case of unsteady motion at either subsonic or super- 
sonic flight speed, the basic equation may be taken in the 


form 
Vir Ya Yw Yr =0 (18) 


where t=a;,t’, do is the speed of sound in the undisturbed 
region of the field, ¢’ is time, and wis the perturbation velocity 
potential or any of the perturbation velocity components. 
Equation (18) is the acoustic equation for small disturbances 
in three space dimensions and holds for a system of Cartesian 
coordinates fixed relative to the undisturbed air. In appli- 
cations to wing theory, therefore, the wings move relative to 
fixed axes. 

In the derivation of a useful theorem it is convenient to 
treat thin wings at small angles of attack and to assume that 
the motion takes place in the zy plane. The visualization of 
the time and geometry relations is relatively easy for two- 
dimensional wings moving at a uniform speed, as indicated in 
figure 2. The airfoil starts at time ¢’=0 and moves to the left 
at a constant velocity Us so that the trace of the leading edge 
in the zt plane is z= — Ud’ = — Mot, and the trailing-edge trace 
is r=2a— Mt. The lines z= +t and s=2a-bi are the traces 
of the extremeties of the regions affected by the acoustic 
waves set in motion at t’=0 by the leading and trailing edges. 
In figure 2, the wing has traveled a time T’=T/a) and the 
boundary condition determining the wing shape during the 
motion will be fixed by prescribing the value of vertical 
induced velocity w over the region “swept out” of the zt 
plane by the wing. In order to determine a reciprocal theo- 
rem, a second wing is assumed to start at the final position of 
the first wing and to move with negative velocity until it has 
réached the initial position of the first wing. With theso 
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7 
FIGURE 2.—T Two-dimensional wing in zt plane. 


concepts in mind, it follows from equation (4) that the 


relation 
D?’ _ fae’ 
JJescSidedt| e 


holds where y and y’ are, respectively, the perturbation po- 
tentials for the forward and reverse motions. The region of 
integration is determined by the area occupied in the zt 
plane by the wing, and the Kutta condition is assumed to 
apply to the trailing edge when in subsonic flight. 

If the left side of equation (19) is integrated by parts, the 
general relation becomes 


Op, aay Og" l 


If the motion of a three-dimensional wing is to be studied, 
equations (19) and (20) must be modified to include an inte- 
gration with respect to y. 

Two further changes in equation (20) serve to simplify 
applications. In the first place, asymmetry is restored to 
the expression if two distinct systems of axes are used as in 
equation (14); in the second place, the pressure relation 


Ò Ò 
P- P= — posn=— Polo; 


(19) 


(20) 


(21) 


where p denotes undisturbed pressure, permits the introduc- 
tion of pressure p in the integrands. The final expression 
for the three-dimensional case is, therefore, 


ff Pr (a1, Yr, by) Wa (21, Yi, td dxidy dt, 


= |f [onun tolen vntddzadyadt, (22) 


where the two motions now follow the same path in reverse 
directions but are referred to the two sets of oppositely 
oricnted axes satisfying the relations 


{j= — r-H E, Y= —yY:+n, 21) hi=— htr (23) 


where ¢, 7, and 7 fix the relative positions of the two origins. 
Figure 3 indicates one possible orientation of the axes. 
Equation (22) reduces to a much simpler form, provided 
further restrictions are imposed on the upwash functions 
321605—b5——18 
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FIGURE 3.—Coordinate systems for three-dimaensional wing in unsteady motion. 


We (21, Y1, h) and w(t, Ye, &). In order to fix the idea, con- 
sider the case in which the wings have traveled a time T”, 
(T'=T/a) and a distance U)7”’=M,T. Let the two sys- 
tems of coordinate axes be placed such that t;=0 sets the 
starting time of the forward motion and t=0 sets the start- 
ing time of the reverse motion; the two origins are further- 
more oriented such that they are at opposite ends of the 
root chord of the common plan form. Equation (22) then 
becomes 


T 
f at, f Í Pir Yis tD Wa Bi, Yi, t)dx,dyy 


Pı (t) oe 
=|, zf l Pa (a, Ya, ta); (£2, Yo, te) dead ys 
/ A ta) 


where the functions w, and w have an implicit dependence 
upon T. If w, and w remain constant for z,+Mpf,—=const. 
or %2+M)f,=const., the expression 


fp (a, Yi, T) ws (21; Yis T)dæzdy 
P(T) 


= f | rales Ya, L)Wi(L2, Y2, F)dtady2 (24) 
P(T) 


follows after taking a derivative with respect to T of the 
original equality. Equations (24) and (14) are now equiva- 
lent in form, with T taking the role of an auxiliary param- 
eter. In this way, certain classes of unsteady motions can be 
treated simultaneously with steady motions. — 

In the applications to follow it will be convenient to intro- 
duce into equation (24) upwash functions of the indicial 
type; that is, functions that are zero up to a fixed time and, 
after experiencing a finite discontinuity, remain constant 
for all subsequent values of time. Such indicial or step 
variations can be assumed, say, for angle of attack, rate of 
pitch, and rate of roll since they satisfy the requirements 
underlying the derivation of equation (24). This choice of 
functions will prove to be advantageous in that the integrals 
of the responsive pressures will yield results relating the 
wing characteristics. Theorems to be given later will speak 
specifically of steady and indicial motions. Itis to be under- 
stood, however, that the indicial results can be further 
extended when the same wing is assumed to be executing 
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identical motions in forward and reverse flight. Thus, by 
means of Duhamel’s integral (see, e. g., ref. 20), if f(t) is 
the response in the wing characteristic to a step variation 
in w at time t=0, the response to an arbitrary variation 
with time of w can be written 


FO= 5 [Fe—Dw(r)dr (25) 


If it is known, for example, that the lift per unit angle of 
attack is the same at corresponding values of time for a 
wing experiencing an indicial angle-of-attack change in 
forward and reverse flight, it follows that the build-up of 
lift is the same at corresponding values of time for all for- 
ward and reverse motions, provided the time histories of 
the motions are the same. The equivalence of lift would 
thus be established, for instance, for oscillatory variations 
in angle of attack. 

An alternative study of reverse-flow theorems for oscilla- 
tory motions could be based upon the modified wave equation 


Xat Xw t Xato x= 


which results from setting. y=e''x in equation (18). Such a 
study would corroborate the conclusions drawn from equa- 
tions (24) and (25). 


APPLICATIONS 


The results of the foregoing analysis may be employed to 
determine a number of special theorems that are particularly 
useful, in the calculation of the aerodynamic characteristics 
of twisted wings and of wings in nonuniform downwash fields. 
The theorems apply equally to wings acting either alone or, 
In certain cases, in combination with other wings or with 
cylindrical bodies having their generators alined with the 2 
axis. Moreover, they apply not only to wings in steady 
motion but also to wings performing unsteady motions of the 
indicial type, or unsteady motions derivable therefrom. For 
wings in more complex unsteady motions, however, it will 
be necessary to refer to the more general results of equation 
(22). Some problems of this nature will be described at the 
end of this section. 

The applications to be included are exact within the frame- 
work of linear theory and involve no further rectrictions on 
the wing plan form or Mach number except in certain indi- 
cated cases where it will be convenient to use results based 
on slender-wing theory. The examples are intended to be 
representative in nature. 


REVERSAL THEOREMS—STEADY AND INDICIAL MOTIONS 


Reversal theorems are defined here as relations between 
the aerodynamic characteristics of identical wings executing 
the same type of motions in forward and reverse flight. The 
results presented in this section apply not only to single 
wings in steady motion but also to combinations of wings, as 
in cascades or multiplanes, performing either steady motion 
or motions of the indicial type. 


DRAG OF SYMMETRICAL NONLIFTING WINGS 


The drag of a symmetrical, sharp-edged wing in linear 
theory may be determined by integrating over the plan form 
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the product of the pressure and the slope in the z direction 
of the wing surface; when the wing has blunt edges with 
slopes having square-root singularities, these singularities 
yield an added contribution (ref. 19). In general, therefore, 
the drag D of a symmetrical section is given by 


D=D,+2 iH} ; (F ds (26) 


where D, is the drag attributable to the edges. 

If the subscripts 1 and 2 refer to the same wing in forward 
and reverse flow, respectively, and with the two systems of 
axes introduced in equation (23), local slopes are related as 


follows 
dz; (21,Y1,L’) 


/ 1 


_ dala T) 
= an (27) 


u 


Equations (26) and (27), together with the reciprocal relation 


(24), yield 
D—Di=f Ja ns ($2) as.=— i p: ($2) asi= 


ffa fa) 


Since the geometry of the wing is fixed, the edge contribu- 
tions are the same, 
(De) = (De)s (29) 
and, consequently, 
D, =D, (30) 
which confirms the relation stated in reference 9. 


THEOREM: The pressure drag in steady or indicial 
motion of symmetrical nonlifting wings is the same in 
forward and reverse flight. 


LIFT ON FLAT-PLATE WINGS 


The lift Z of a wing may be determined by integrating tho 
differential pressure Ap=Pı— Puy ‘over the wing plan form, 


thus 
Te f { ApdS 
P 


For flat-plate wings, the local angle of attack of the wing sur- 
face is a constant 


(31) 


Mg (x2, Yzy T) =a, = const, (32) 
Application of equations (31), and (32), and (24) yee the 


following: 
a nr hans = La 


Dy, /o,=L4/ 02 (33) 


THEOREM: The lift per unit angle of attack of flat- 
plate wings in steady or indicial motion is the same in 
forward and reverse flight. 


Cy (x4, Yis T) =a,=const., 


or 
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This theorem generalizes the relation previously given by 
Brown (ref. 5) for steady motion. 


DAMPING IN ROLL OF FLAT-PLATE WINGS 


The rolling moment L’ exerted on a wing, following the 
usual sign convention, is given by 


Das f f yApdS 
P : 


The local angle of attack due to rotation about the z axis is 


(34) 


M'Y Dal Ya piy t 
i Uo : i Uo Uo (35) 


where p’ is the angular velocity of roll, assumed constant. 
Application of equations (34), (85), and (24) yields the 
following: 


i , / ` 
P 1 F 1 
pi Ly 


= BL —_ Pr _ 
Bf yApdS2 Ti, | ysdpsd8, Ti 
f; P} 


or 


Ly [pi =L |p’ 


THEOREM: The rolling moment per unit angular 
rolling velocity of flat-plate wings in steady or indicial 
motion is the same in forward and reverse flight. 


(36) 


DAMPING IN PITCH OF FLAT-PLATE WINGS 

Consider a wing, first, in forward flight and pitching with 
a uniform angular velocity qı about a lateral axis; second, 
in reverse flight and pitching with angular velocity q: about 
another lateral axis. Place each wing in a coordinate system 
such that the y axis coincides with the axis of rotation and 
designate the distances to the moment axes by 2 with 
proper subscripts, as shown in figure 4. Jn such a coordinate 
system, the pitching moment M- exerted on a wing, following 
the usual sign convention, is 


ae f f E (37) 
P 


FIGURE 4,—Coordinate systems and symbols used in discussion of reverse-flow theorems 
for pitching wings. 


The local angle-of-attack distributions due to rotations 
about the y axis are 
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a= Fp Ti; a= FF a (E—z) (38) 
Application of equations (37), (38), and (24) yields the fol- 
lowing: 


iP (Ma FE— tala =f i (ta — rH — tap dS = 


IA f -zand = | f T:4pıd S1 = 
P, Pi l 


g Jy mA pds =E KMar + E— 200) La) 
2 


or 
(M 1)zo1 teca M 2)z + (E—202)L2 (39) 
1 q2 A 


This equation indicates that the pitching moment due to 
pitching velocity is, in general, not the same for wings in 
forward and reverse flight. However, if =ze=—é, the 
pitching moment per unit angular pitching velocity of flat- 
plate wings in steady or indicial motion is invariant. 


SPECIAL RECIPROCAL THEOREMS AND APPLICATIONS 


In the following section, severul special reciprocal theorems 
will be derived and applications will be illustrated. Recipro- 
cal theorems, in contrast to reversal theorems treated in the 
preceding section, are defined here as relations between the 
fxerodynamic properties of wings in forward and reverse flight 
that have dissimilar camber, twist, and thickness distribu- 
tions but have the same plan forms. The motions may or 
may not be similar, although it is assumed in this section 
that both wings are in either steady motion or unsteady 
motion of the indicial type. As noted in the preceding sec- 
tion, the results apply equally to wings acting alone or in 
combination. 


‘SYMMETRIC NONLIFTING WINGS—STEADY MOTION 


The problems of paramount interest in the application of 
the general relations are found from considerations of pres- 
sure integrals over. lifting surface; such problems will be 
given detailed treatment later. In the present section, a 
brief indication is given of the manner in which useful results 
can be derived for symmetric wings. The discussion will 
be limited to steady-state, two-dimensional, subsonic pres- 


sure fields although fairly obvious extensions can be carried 


out. 

If the geometry of a real symmetric airfoil is prescribed 
the theoretical pressure distribution exists and is unique. 
Lf, however, the pressure distribution is prescribed, a real 
airfoil does not necessarily exist, but by means of reciprocal 
relations it is possible to derive certain conditions of com- 
patibility that need to be imposed. Consider, therefore, 
the two subsonic solutions ~ 


wlz) =U Jr 5; —a<n<a 


Un (ta), 


u(z)=0, 


and 
Wa (£a) 


ce a ee ee ee me a a 
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The first solution has square-root singularities in w at each 
end of the airfoil-and, correspondingly, singularities in u occur 
just ahead of the point s=—a and just behind z=a. On the 
other hand, equation (14) certainly applies if we, is zero at 
z=-ta. If the origins of the two systems of axes are at 
the same position, it follows from equation (14) that w(x) 
must satisfy the relation 


fl 


This result is useful in the calculation of airfoil shapes in- 
volving a change in pressure distribution from that of a known 
reference profile. ‘The restriction on w at the nose and tail 
implies that the derived and reference profiles have the same 
slope and radius of curvature at those points. The restric- 


(40) 


tion on t, as given in equation (40), can be interpreted as a | 


condition that must exist by virtue of the fact that the drag 
of an airfoil in two-dimensional potential flow in zero. 
As a second example, consider the solutions 


tU, l — trt Uo . 
ne =a ned =; —a<n< 


that represent velocity and slope of a thin ellipse of .thickness 
t and chord 2a. If wis chosen as above, such that it van- 
ishes at the nose and tail of the airfoil, if u is the correspond- 
ing velocity distribution, and if the two sets of axes are as 
before, equation (14) yields 


f: dza gp —— (° Ua@a)taddes 
ada ° -a ya — t 
From tbis result, together with the general closure condition, 


(41) 


a necessary condition for the closure of the second airfoil is 


* Up (ta)t2dry o 


ae (42) 


—a 


As a final example, consider the solutions for —a<2,<a 


“9 


2tUo 
tı) = a—2r 
uak 1) 343a? ( 1); 
Sm me — n2 
W, (2) = 2tUo 421 atı — a 


3 V3 a? y a?7— x7 i 


representing velocity and slope of a thin Joukowsky type 
airfoil. In this case, w, vanishes at the tail and the down- 
wash distribution w for the reverse wing may have & square- 
root singularity at the nose. The nose of the first wing is, 
however, blunt and for equation (14) to apply the second 
wing must have a cusped tail. Under these conditions, 
equation (14) yields 


2tU,? f° dz: 

33a? e. 221) dat ati 
__ 2tU, (° 22°—ax,—a 
-2 f a ulada 


ped 


| 
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Making the substitution 2,;—=—2, in this equation and inte- 
erating the left side by parts, one has 


To] enaz 2f" zazidan? 


* 


¢ Q27,7+-a2,—a? 
a eens a A 
a gan Uata) dx. 


For all real airfoils with cusped trailing edges, therefore, the 
area, A, can be expressed as 


a 2297-+ a%y—Q? Ua(2a) des 
~a a—z,? Uo 
LIFT—STEADY AND INDICIAL MOTION 


The reciprocal theorems offer considerable advantage in 
the calculation of the lift of wings having & nonuniform 
angle-of-attack distribution or of wings in a stream having 
nonuniform flow directions. For these applications, it is 
convenient to consider a special form of the reciprocal theo- 
rem which relates the lift on a wing having arbitrary distri- 
bution of local angle of attack to that of the flat-plate wing 
of identical plan form in flight in the reverse direction. 
Since the solution of this latter problem is often known or 


A= Kg (4 3) 


‘can be found relatively easily, the solution of the original 


problem is facilitated in' many instances. 

Lift of arbitrarily cambered wings.—Consider two wings of 
identical plan form in flight in opposite directions, as shown 
in figure 5. Wing 1 is arbitrarily cambered and twisted and 


wing 2 is flat. 
(44) 


Oy = (xy Yis T) » Q3 const. 


FIsurRe 5.—Coordinate systems and symbols used in discussion of rolation betweon lift of 
arbitrarily cambered wings and loading on flat-plato wings. 


Application of equations (44) and (24) yields the following: 


æla = ff aanpas 1 = fesdpedSs 
P, P 


2 
L= f fo (22: dSa 
P: 


THEOREM: The lift in steady or indicial motion of a 
wing having arbitrary twist and camber is equal to the 
integral over the plan form of the product of the local angle 
of attack and the loading per unit angle of attack at the 
corresponding point of a flat-plate wing of identical plan 
form in flight in the reverse direction. 


(45) 


RECIPROCITY RELATIONS IN AHRODYNAMICS 


Equation (45) may be used to derive Munk’s integral for- 
mula for the lift of an arbitrarily curved airfoil in subsonic 
flow. Consider airfoils 1 and 2 placed in their respective 
coordinate systems, as indicated in figure 6. The angle-of- 
attack distributions on the two airfoils are given by ` 


(46) 


FIGURE 6.—S8ketch of arbitrarily cambered airfoil Dlustrating symbols used in equations (46) 
through (50). 


The loading per unit angle of attack on airfoil 2 is 


ate, 


a—z, ` 


Ap: 440 [attt 


a B atx &8 


where qo is free-stream dynamic pressure ($p9U/,?). Substi- 
tution into equation (47) yields the lift formula 


(47) 


L= f" dz, Ze om (48) 


B J-adx,V a—2, 


The corresponding formula for the lift of a tier of curved 
airfoils may also be derived similarly from the expression 
for the loading on an equivalent tier of flat airfoils. Con- 
sider, for example, an unstaggered lattice of flat-plate 
airfoils arranged vertically. If the gap distance between the 
Dlates is h, the loading per unit angle of attack is 


APs 40 soch (32 [sinh [aB] 

a B NBh) V sinh [x(a+23)/Bh) 
The formula for the lift on one of a lattice of identically 
cambered airfoils ? is therefore 


___ 490 (=) f dz, /sinh [r(@+21)/Bh] 
L= ch (Bh) J-a da V sinh fala — zih] - 
The load distribution per unit angle of attack for a two- 


dimensional supersonic wing is 


APs _ 44 


ao 8 (49) 
and, from equation (47), the lift is 
— 290 (* dz Í 
EH B f dz, an (50) 


The extension of this result to include supersonic-edged wings 
straight trailing edges leads to a result given originally by 
Lagerstrom and Van Dyke (ref. 21). If, as in figure 7, the 

3 This result, as well as the detailed pressure distribution, has been derived by Mr. Paul 
F. Byrd of the Ames Acronautical Laboratory by means of a direct inversion of the singular 


integral equation relating the aerodynamic pressure and geometry of the airfofl. His work 
has also been used in deriving equations (80) and (81). 
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FIGURE 7.—SkKetch of supersonic-edged wing Mlustrating symbols used in equations (51) 
and (&2). 


sweep angle of the straightedge is A, the load distribution 
per unit angle of attack of the reversed wing is 


AD £ qo cos A 


Se 51 

a W1—M,? cos? A Po 

and the lift is | 
— 4go dz, 

L= (52) 


-l —— dzd 
sec ?A— Mè 5 da, h 
1 


À less obvious application yields the build-up of lift with 
time of an arbitrarily cambered two-dimensional supersonic 
airfoil starting impulsively from rest at a constant speed. 
Rewriting equation (45), lift is 


= s-M0T da (xa, T) Ap: (£T) dite 


eee 


—MoT dx, a3 


(53) 


Figure 8 presents sketches showing the final positions of the 


fi 


2a 2 a- 2a 2a ~ 
ISEE TEM tM l < Mel Mi T 


FIGURE 8&.—Regions in which different forms of equation (54) apply. 


airfoils relative to each other for various values of T. In 
reference 22 the expressions for the loading in equation (53) . 


_ are given. Over the intervals denoted by ©, ©, and ©) in 


the sketch, these expressions are 


; A 4 
Region () aa A 

e ADs __ Ado Moz +T 
Region © a AM are cos ZI AAT m 

hails Pan A : st arc sing) | 

; A 4 

Region @ 2 = T 
= 


Lift on a wing in a nonuniform downwash field.—The 
reciprocal theorem of equation (45) can also provide a par- 
ticularly good method of determining the lift on a wing in 
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certain nonuniform downwash fields of known structure. 
Such problems arise whenever a wing acts in the presence 
of other wings, bodies, or propellers but is always of prime 
concern in the determination of the lift on a tail acting in 
the downwash field of a wing. In most problems, the down- 
wash velocities at the position of the tail may be considered 
to be constant in the longitudinal direction and to vary in 
the spanwise direction, thus 


kaji 


a= Yn), a=const. (55) 


ref fa (2) dem f(a 
P, 


where l is the span load distribution associated with the 
load distribution Ap:. Summarizing, the lift in steady or 
indicial motion of a wing in a downwash field which varies 
across the span is equal to the integral over the span of the 
product of the local angle of attack and the span loading per 
unit angle of attack at the corresponding spanwise station 
of a flat-plate wing of identical plan form in flight in the 
reverse direction. This statement generalizes the result 
given recently by Alden and Schindel (ref. 23) for steady 
flow about wings having supersonic leading and trailing 
edges and streamwise side edges. 

As for example, consider the problem of determining the 
lift on a wing at a geometrical angle of attack of zero resulting 
from the presence of an infinite line vortex of strength T 
extending in the flight direction. The wing will be considered 
to have such a plan form that its span loading when in 
flight in the reverse direction is elliptic. The notation is as 
shown in figure 9. For this problem, therefore, the span 


and 


(56) 


Fiaure 9.—View of wing and neighboring vortex. 
loading of the wing in reverse flight is given by 
l = 214 1 yr 

-_ S 


Q2 TAEA 


(57) 


The local angle of attack of the original wing due to the 
presence of the vortex is given by 


T ` YB 


o Aah i — BP 129) 
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Substitution of equations (57) and (58) into equation (56) 
yields the following formula for ‘the lift: 


By, apr 
a) tt 


F- 


The lift on a wing in the vicinity of a number of such 
vortices may be found by superposition. 

Lift due to deflection of a portion of the wing surface.— 
Let a portion P’ of the surface of wing 1 be deflected a con- 
stant angle ô and the remainder of the wing, be a flat plate 
alined with the free-stream direction.’ Let wing 2 be a 
flat-plate wing inclined at an angle of attack a, thus 


ô on P’ 
Q = 


Qg = const. (6 0) 


0 elsewhere 
Substitution of equation (60) into (45) yields the following 


result: 
L f ((Aps 
am Ca) 
2 


The lift in steady or indicial motion per unit angular 
deflection. of a portion of the wing surface is thus equal to 
the lift per unit angle of attack on the corresponding portion 
of a flat-plate wing in flight in the reverse direction. This 
generalizes a result given previously by Morikawa and 
Puckett (ref. 24) for steady flow about low-aspect-ratio wings. 


This rule is very useful in the determination of the lift 
resulting from the deflection of a flap or control surface. 


(61) 


‘This is particularly true for supersonic speeds since the load- 


ing on the related flat wing is often a constant over a largo 
portion of the area of integration. : 


As a further. example, consider the case of a low-aspect- 


‘ratio wing having a straight trailing edge and mounted on 


an infinite cylindrical body of revolution. The entire wing- 
body combination is at zero angle of attack except for the 
flaps on the rear of the wing that are deflected an angle ô. The 
problem is to determine the lift on the entire wing-body 
combination due to the deflection of the flaps. Slender- 
wing-theory results of reference 25 are to be used. The 
notation is indicated in figure 10. 


x 


Fiavre 10.—View of lifting slender wing-body combinations. 


RECIPROCITY RELATIONS IN AERODYNAMICS 


The solution of this problem is particularly facilitated by 
the fact that slender-wing theory indicates that the loading 
on wing 2 is concentrated on the leading edge, as shown in 
the figure. Therefore, the lift of wing 1 is found by inte- 
grating the span loading curve of wing 2 over the portion of 
the span between y=f and y=s. Thus i 


e- SSE Ea 
flap 
The span loading on wing 2 is given by 


sont (0) 


The lift due to the deflection of the flap is therefore 


2 4 
Lop =4q080") (1-5 a-y- (5-4 i 


1 at . 1—2(f/ + (afs)! 
aO Ha) are sin = 


a? -B +(a/9)9 (s) —2(a/8)! 
see FSG 


A plot of the results is shown in figure 11. The lift per unit 
angle of flap deflection (Z/6),1., has been nondimensional- 


> 


(62) 


(63) 


(64) 


(278) iop/(L/a)w-8 


a/s 
FIGURE 11.—Lift of a slender wing-body combination resulting from flap deflection, 


ized by dividing by the lift per unit angle of attack (L/a)w-s 
of a slender wing-body combination of identical plan form. 
From reference 25, (L/æ)w-»s is given by 


(L/e) 7 p= 2 gos? (1 —5) | 


s 


(65) 


ROLLING MOMENT—STEADY AND INDICIAL MOTION 


The calculation of the rolling-moment characteristics of 
wings having a nonuniform angle-of-attack distribution or 
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of wings in a stream with nonuniform-flow directions can be 
performed in many cases through use of the reciprocal 
theorem in a manner analogous to that described for the 
lift characteristics in the preceding section. In every case, 
the rolling moment of the given wing will be related to the 
lift on a rolling flat-plate wing. For the sake of simplicity, 
all the present examples will be confined to the case where 
the rolling moments are evaluated about the x axis, con- 
sidered to lie in the plane of symmetry. 

Rolling moment of arbitrarily cambered wings.—Consider 
two wings of identical plan form in flight in opposite direc- 
tions, as shown in figure 12. If the local angle-of-attack 


FIGURE 12.—Coordinate systems and symbols used in discussion of relation between rolling 
moment of arbitrarily cambered wings and the loading on rolling flat-plate wings. 


distribution of wing 1 is arbitrary and that of wing 2 varies 
linearly with y (which might be likened to either a wing with 
linear twist distribution or to a flat-plate wing rolling about 
the z axis with constant angular velocity p,’), 


- 7 
a, =a (tY T), te (66) 
0 


and the following relations can be written: 


tim JF napas | | napas = Sy a a2,Ap dS = 
; l l 1 1 
Uo = APs 
J | jf tedS.= JJ a (Fir) dS: (67) 
2 2 


THEOREM.: The rolling moment in steady or indicial 
motion of a wing having arbitrary twist and camber is 
equal to the integral over the plan form of the product of 
the local angle of attack and the loading per unit (p,’/U)) 
at the corresponding point of a rolling flat-plate wing of 
identical plan form in flight in the reverse direction. 


Applications of this theorem follow in a manner very 
similar to that described previously for the- corresponding 
theorem regarding lift. 

Rolling moment on a wing in a nonuniform downwash 
field.—Consider a wing placed in a flow field in which the 
downwash velocities at the position of the wing are constant 
in the longitudinal direction and vary in the spanwise direc- 
tion. The related wing is again a flat-plate wing rolling with 
an angular velocity p’ as described in the preceding section, 
thus 


Oy = Oy UT), a (68) 
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The rolling moment of the first wing is then given by 


[fe Gin) fi a Gi) o 


or, in words, the rolling moment in steady or indicial motion 
of a wing in a downwash field which varies across the span is 
equal to the integral over the span of the product of the local 
angle of attack and the span loading per unit (p.’/U ) at the 
corresponding spanwise station of a rolling flat-plate wing of 
identical plan form in flight in the reverse direction. 

Rolling moment due to deflection of a portion of. the wing 
surface.—Let a portion P’ of the surface of a wing be de- 
flected a constant angle 6 and the remainder of the wing be a 
flat-plate alined with the free-stream direction. The related 
wing is a flat-plate wing rolling with angular velocity p,’ 

f 
5 -| on P ay Pe Ys 
Uo 


(70) 
0 elsewhere 


Substitution from equation (70) into (67) yields the following 


result: 
3G) om 


Thus, the rolling moment in steady or indicial motion due 
to a given angular deflection of a portion of the wing surface is 
equal to the lift per unit (p2’/U,) on the corresponding portion 
of a rolling flat-plate wing of identical plan form in flight in 
the reverse direction. 

As an example, consider a wing-body combination con- 
sisting of a low-aspect-ratio wing having a straight trailing 


edge mounted on an infinite circular cylinder, as shown in . 


figure 13. The body is at zero angle of attack, the right wing 


Uo 


FIGURE 13.—S8ketch illustrating symbols used in discussion of rolling moment resulting 
from differential deflection of wings of slender wing-body combination. 


P” is deflected an angle ô, and the left wing P’’ is deflected —é. 
The problem is to determine by means of slender-wing theory 
the rollmg moment exerted on the entire wing-body com- 
bination. The notation is indicated in figure 13. 


(72) 


0 on body 
+ô on right wing, P’ 


R on left wing, P” P 
__ P2 yz 

ay = Q = 

Uo 


Since slender-wing theory indicates that the loading on wing 


2 is concentrated on the leading edge, the rolling moment of 
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oe 1 can be found by integrating the span loading on wing 


z = = (I0 26.) asa J" Grig,) at 
perm a f, GAT) 2 


where, from reference 26, 


hb 
pa'] Uo 


2 2a8 a a! ; as 
=a (1 T arc cos Fea (yt =) vay — at 


Z =(y -2Y arc cosh EE 
(73) 
The resulting expression for the rolling moment is 
Ly =4qoôs* Fío, k) Bee 2 m+) +2 arc COS 7755 Sop = 
ak “(1 -P |+E (y, k) È (G- 2+ (1 + 


2 oR \, 2 
= are cos zry) += RU si LE In 5 et 


R 2 2R 
= (1—2) (1+2 are, cos = z) (74) 
where - 
R=a/s 
g= aoin 
k=y/I 


A plot of the results is shown in figure 14. The rolling 
moment has been nondimensionalized by dividing by tho 
value corresponding to that of the wing alone (R=0). 


PITCHING MOMENT—STEADY AND INDICIAL MOTION 


A number of useful relations regarding the pitching- 
moment characteristics of wings may be found through 
application of the reciprocal theorem. Since the general 


‘procedure is closely analogous to that of the preceding sec- 


tions, the following discussion will be brief. 

Pitching moment of arbitrarily cambered wing.—Consider 
the problem of determining the pitching moment Mı about 
the origin of wing 1 possessing an arbitrary distribution of 
camber. The related wing in flight in the reverse direction, 
wing 2, is a flat-plate wing of identical plan form pitching 
about the moment axis of wing 1, as indicated in figure 15, 
thus 


az= = (75) 


eg\=a(,Yil), ni- £) 


RECIPROCITY RELATIONS IN AERODYNAMICS 


o/s 


FIGURE 14.—Rolling moment of slender wing-body combination resulting from differential 
deflection of wings. 


FIGURE 15,—Coordinate systems and symbols used in discussion of relation between pitching 
moment on arbitrarily cambered wings and the loading on pitching flat-plate wings. 


The pitching moment of wing 1 is given by 


mom [fonts [fa (as 
P 1 P 1 
= Apa 
~ J ü G Tm) eo) 
2 


THEOREM: The pitching moment in steady or indicial 
motion of a wing having arbitrary twist and camber is 
equal to the integral over the plan form of the product 
of the local angle of attack and the loading per unit 
(q2/U>) at the corresponding point of a flat-plate wing 
of identical plan form in flight in the reverse direction 
and pitching about the moment axis of the first wing. 
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The necessity for pitching wing 2 about the moment axis 


‘of wing 1 may be removed by considering wing 2 to be 


re-expressed in terms of two component wings having angle- 
of-attack distributions given by 


da= qatay — Tos’) a Ulta £) 
Uo — To 


Wing 2’ is thus pitching with angular velocity q} about an 
axis at %2=2y and wing 2” is a flat-plate wing at a con- 
stant angle of attack. The pitching moment on wing 1 is 


then given by 
Apa 
Sa | fo Ge 
Pon 


i Apr 
oe ffa Cr; Uo 
Py 


= f fou (F7) dS tew | f oa (AP) dS» (78) 
Po: Py. f 


Applications of pitching-moment theorem.—The applica- 
tion of equation (76) or (78) to problems analogous to those 
discussed in the preceding section can be carried out in a 
straightforward manner. Consider, first, unstaggered lat- 
tices on airfoils such that the airfoils in lattice 1 have arbi- 
trary camber distributions and those in lattice 2 are flat 


=const. (77) 


aS gr 


plates pitching about their midchord positions. The angles 
of attack in the two lattice systems are 
d 
on(21,21)= ae 2 (%2,22)= T > (79) 


and the load distribution on each airfoil in lattice 2 is, in 
subsonic steady flow, 


Ap(£)= Go (=) arc cos | sech (=) cosh (FR) | (80) 


where 2a is chord length. Equation (76) yields, for pitch- 
ing moment of the first airfoil about its midchord pomt, 
the result 


M — 41A f Gar iig cos| sech (=) cosh (5) | dz, (81) 


K —a@ ax 


A second example, illustrating unsteady effects, is the fol- 
lowing: Let wing 1 be a flat-plate wing, then a; is constant, 
and equation (78) simplifies to 


Mı La: 


M E 
a1 =z T e £) 


0+ 


(82) 


where J,’ is the lift on wing 2’ pitching about ty? =z", and 
L,’’ is the lift on an inclined flat-plate wing. Equation (82) 
may be expressed in terms of conventional stability deriva- 
tives as follows: 


(On = (Court (72 —* Or) (83) 
An application of this result to unsteady-flow problems is 
indicated in figure 16 obtained from indicial-lift and pitching- 
moment results of reference 27. This figure shows the 
growth of lift and pitching moment on triangular wings with 


Sa ey 
ee ee ee „n a terre 


$ 15 zo 
Chord lengths traveled 


Fiaure 16.—Indiclal lift and pitching moment on triangular wings with supersonic edges. 


supersonic edges at a Mach number of 2 following indicial 
angle-of-attack and pitching-velocity changes. In these 
results, the rotation and moment axes are always at the 
leading edge or apex, therefore, t:’=0 and =c. It may 


be seen that the three curves are related in the simple linear © 


manner indicated by equation (83). 

If a, is independent of z and varies only in the spanwise 
direction, that is, if a=«,(y,), the pitching moment on 
wing 1 is given by the following equation, analogous to 
equation (56) for lift: 


Tf a portion P’ of the surface of wing 1 is deflected a con- 


stant angle 5 and the remainder of the wing is a flat plate 


alined with the free-stream direction, the following relations 
hold: 


p on P’ . 
0 elsewhere 
and 


ff ( Ap 
ô ga'[U 9 
P'y 


Sytari) || (FB) dSx (88) 
Pe > 
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RECIPROCAL RELATIONS INVOLVING MOTION INTO A GUST 


All previous applications that have been considered were 
derived from equation (24). In the present section, the 
more general equation (22) will be used to develop two the- 
orems which relate the build-up of lift on a wing entering 
a gust and the build-up of circulation on the same wing 
moving indicially but in the opposite direction. The rela- 
tions to be obtained hold for the Mach number range for 


which the wave equation applies. Under the special assump- 


tions of incompressible flow, the results in two dimensions 
establish a direct connection between the circulation function 
calculated by Wagner (ref. 28) and the gust lift curve cal- 
culated by Kiissner (ref. 29) and von Kármán and Sears 
(ref. 30). A proof of the connection between these functions 
for two-dimensional incompressible flow has been given by 
Sears (ref. 31). 


TWO-DIMENSIONAL FLOW 


A flat plate is assumed to be moving in two modes of 
motion: In the motion associated with the 2:,2,f,, axes, the 
wing starts at time zero (4=0) and moves at a constant 
velocity Uo and at a constant angle of attack; the motion 
associated with the 2,,2,,t, axes starts at time zero (¢,=0) 
with the wing moving‘in the opposite direction at a velocity 
U, and entering a sharp-edged gust. The gust exists for 
all x, less than zero and has a vertical velocity w=—a, Uo. 
The two wings, therefore, have angles of attack as follows: 


a(z, t) =a =const. for—Myk<m<2a—Moh, 4>0 
—Moti<a<0 0 s é s 2a/Mo 
Cty (xy f t) =a,—=const. fo 
: — Moti Lt: <L 2a— Mot 2a] Mo sS éi 


The two-dimensional form of equation (22) yields 
T — Mot > 
aa fds S APdEn thde1= potoa | | a 2nd ag 
A (86) 


where the region A is bounded by the lines z,=0, xı=2a— 
Mt, L= T, and zı=— Md. The integral on the right can be 
rewritten as a line integral by means of the identity 


— Í P cos t, a)ds= | E as 


and equation (86) becomes 


a f E T f tace et 
g 


4 


where the line integral extends around the boundary of the 
region A. Since Ag; vanishes on the lines tz= — Mog, and 
t=0, the equation becomes 


fr (tdi, = a A (2 20 ayy 
Qe ae 1 Podots 2a—MoT Y2\ La, M, 2 


Differentiation with respect to T yields 


onL,(T) =poUoaAge(2a—Mof, T) 


RECIPROCITY RELATIONS IN ABRODYNAMICS 


The discontinuity i In gis evaluated at the trailing edge at 
time T and is therefore equal to the circulation I, of the air- 
foil at time T. The equality thus becomes 


aD — p72) (87) 


THEOREM: The circulation per unit angle of attack 
of a flat plate moving indicially with a velocity Uy is 
proportional to the lift per unit a, of the plate entering 
a sharp-edged gust having a uniform vertical velocity 
equal to wg=—a,Up. 

In figure 17, the time variation of these variables, as well 
as the lift of the indicial wing, is indicated for low speed and 
for flight Mach numbers equal to 0.8 and 1.46 as determined 
from references 22 and 27. 


5 lo 
Half-chords traveled 
Mo*0.8 ` 


50 5 10 15 
Half-chords traveled 
Mp =1.46 


5 
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FIGURE 17.—Growth of ¢;, and cr with chord lengths traveled. 
C 


THREE-DIMENSIONAL FLOW 


The extension of the above results to three dimensions 
follows directly. The origin of the x, Yı, 21, tr axes is assumed 
to be initially at the foremost point of the wing in reverse 
motion. The two wings have, respectively, angles of attack 
«,=const. over the reverse moving plan form for all values 
of time and o,—const. over the region occupied simultane- 


ously by the forward moving wing and the gust. Equation 
(22) gives 
T 
aa | at, ff Ap,(21, Yi; hjdzidyı = 
P(t) 
petra, | | f a Senni Bedydd (88) 


The integral on the right can be rewritten as a two-dimen- 
sional surface integral by means of the identity 


ff P cos (t, mas—|{ {> dzdydt 


and equation (88) becomes 


T 
o | Llt)dti= podoa, If Ags cos (h, n)dS2 


where the integral on the right extends over the boundary of 
the volume in z2, Y2 f2 space occupied by the wing and the 
gust. The value of Ag; must, of course, vanish on the lead- 
ing edge of the wing and at 4=0. In order to fix the limits 
of integration, suppose the wing is symmetrical about its 
longitudinal axis and let the leading edge of the forward 


267 
wing be given by the equation 
m=f(LY)—Moh or yi=+8(21+Moh) . 
For the reverse wing and its coordinate system, this edge, 
which is now the trailing edge, is 
%=C,—Mots—f(4y2) or Y= +8 —rt:— Moha) 


where c, is the root chord. The reverse-flow integrals of 
equation (88) then become 


T Cy 8 (Cp—Z) 
VAA TE f da f 
o f of 1) 1= Podoads ¢—MoT 2 E 


Agp | za Y2 ea fen) dY: 


where 8(z:) is the local half-span of the wing. Differentia- 
tion with respect to T yields 


Ags | e= MP, yt LS) dYa 
0 . 
(89) 


The discontinuity in ¢ is thus to be integrated spanwise at 
the rearmost point of the indicial wing; this follows from the 
relation 


Ags | ay, TA apley: 


which fixes the vorticity in the wake of the wing once it is 
shed from the trailing edge. 

It remains to mention the nature of the limits +s(A4,7). 
As shown in figure 18, the.span width of the vortex wake at 
the trailing edge is, during the early stages of the motion, 
dependent on the local span width of the wing. The width 
28(MT) of wake is, in fact, equal to the maximum width of 
the portion of the first wing that lies within the gust. After 


s( MoT) 
calf L)= poa | s(MT) 
—8 (410 


HELO, 
T- <t<T 


FIGURE 18,—Sketch illustrating nature of integration limits in equation (89). 
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sufficient time has passed for the vortex wake of the indicial 
wing to.develop its full span width at the trailing edge, 
s(MoT) becomes & or semispan of the wing. From equation 
(89) one may conclude the following: 


THEOREM: The lift per unit æy of a flat-plate wing en- 
tering a sharp-edged gust having a uniform vertical 
velocity equal to w,=—a,U, is proportional, at each 
instant of time, to the spanwise integral at the trailing 
edge of the vorticity shed- by the same wing moving 
indicially in the reverse direction with a velocity Up. 


As a direct example, this theorem has been used to confirm, 
from a knowledge of the indicial solution, the sharp-edged- 
gust lift of the rectangular-plan-form supersonic wing given 
by Miles in reference 32. 


ÅMES ABRONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
Morretr FIELD, Carr., Feb. 19, 1958. 
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